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Abstract. 

Gravitational waves bring about the relative motion of free test masses. 
The detailed knowledge of this motion is important conceptually and practically, 
because the mirrors of laser interferometric detectors of gravitational waves are 
essentially free test masses. There exists an analogy between the motion of free 
masses in the field of a gravitational wave and the motion of free charges in the field 
of an electromagnetic wave. In particular, a gravitational wave drives the masses 
in the plane of the wave-front and also, to a smaller extent, back and forth in the 
direction of the wave's propagation. To describe this motion, we introduce the 
notion of "electric" and "magnetic" components of the gravitational force. This 
analogy is not perfect, but it reflects some important features of the phenomenon. 
Using different methods, we demonstrate the presence and importance of what 
we call the "magnetic" component of motion of free masses. It contributes to 
the variation of distance between a pair of particles. We explicitly derive the 
full response function of a 2-arm laser interferometer to a gravitational wave of 
arbitrary polarization. We give a convenient description of the response function 
in terms of the spin-weighted spherical harmonics. We show that the previously 
ignored "magnetic" component may provide a correction of up to 10%, or so, 
to the usual "electric" component of the response function. The "magnetic" 
contribution must be taken into account in the data analysis, if the parameters 
of the radiating system are not to be mis-estimated. 



1. Introduction 

The large-scale interferometers for the detection of gravitational waves are approaching 
their planned level of sensitivity [Jj. In the coming years, we may be observing the 
rare, but most powerful, sources of gravitational waves (for reviews, see for example 
[2-7]). These initial instruments are scheduled to be upgraded in a few years time. 
The upgraded (advanced) interferometers will be more sensitive than the initial ones. 
They will enable us to determine with high accuracy the physical parameters of various 
astrophysical sources. The knowledge of the detailed and accurate response function 
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of the instrument to the incoming gravitational wave (g.w.) of arbitrary polarization 
becomes an important priority. 

A laser interferometer monitors distances between the central mirror(s) and the 
end-mirrors in the interferometer's arms. In the existing instruments, the length 
I of the interferometer's arm is significantly shorter than the wavelengths A of the 
gravitational waves which the instrument is most sensitive to. The evaluation of the 
change 61 of the distance between two mirrors is usually formulated as 61 /I « h, where 
h is the characteristic amplitude of the incoming wave. This is a correct answer, but it 
is correct only in the lowest, zero-order, approximation in terms of the small parameter 
//A. The next approximation introduces a "magnetic" correction l^l to 51/1, which is 
of the order of h{l/X). This correction depends on the ratio / /A whose numerical value 
is at the level of several percents for the instruments with / — Akm and typical g.w. 
frequencies v — c/A « l(y^Hz. This contribution may be especially important for the 
advanced interferometers operating in the "narrow-band" mode aimed at detecting 
relatively high-frequency quasi-monochromatic g.w. signals. 

In this paper, we start from the motion of free charges in the field of an 
electromagnetic wave, section |21 Then, we turn our attention to the main subject 
of the paper - motion of free masses in the field of a gravitational wave. In section |31 
the positions and motion of free test masses is described in the local inertial reference 
frame associated with one of the masses. This choice of coordinates is the closest 
thing to the global Lorentzian coordinates that are normally used in electrodynamics. 
The identification of the "electric" and "magnetic" components of motion, as well as 
comparison with electrodynamics, are especially transparent in this description. In 
section 0] we exploit the geodesic deviation equation for the derivation of equations of 
motion and identification of the components of the gravitational force. The usually 
written equation, with the curvature tensor in it, is only the zero-order approximation 
in terms of l/X. This approximation is sufficient for the description of the "electric" 
part of the motion, but is insufficient for the description of the "magnetic" part. In 
the next approximation, which is a first order in terms of //A, the geodesic deviation 
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equation includes the derivatives of the curvature tensor, and this approximation is 
fully sufficient for the description of the "magnetic" force and "magnetic" component 
of motion. From these considerations it follows that the component of motion which 
we call, with some reservations, "magnetic" is, in any case, the finite-wavelength 
correction to the usual infinite-wavelength approximation. In the end of this section 
we compare our interpretation with other analogies cited in the literature. In 
section [S] we switch from the positions and trajectories of individual particles to 
the distances between them. The calculation of distances is especially simple in 
the local inertial frame, since the metric tensor, in the required approximation, is 
simply the Minkowski tensor. However, the conclusions about distances, including 
their "magnetic" contributions, do not depend on the choice of coordinates. We show 
that the universal definitions of distance, based on the light travel time and the length 
of spatial geodesies, lead to the same result in the appropriate approximation. In 
section|Hlwe use the derived results for the construction of the response function of a 2- 
arm ground-based interferometer. It is assumed that a gravitational wave of arbitrary 
polarization is coming from arbitrary direction on the sky. We explicitly identify the 
"magnetic" contribution to the response function and demonstrate its importance. 
Then in section \7\ we formulate the response function in terms of the spin-weighted 
spherical harmonics. In section |S1 we give numerical estimates of the "magnetic" 
contribution to the response function of the presently operating instruments, in the 
context of realistic astrophysical sources. We show that the lack of attention to the 
"magnetic" components can lead to a considerable systematic error in the estimation 
of parameters of the g.w. source. Some technical details are described in Appendices. 

2. Motion of free charges in the field of electromagnetic wave 

To understand the origin of the gravitational "magnetic" contribution to the motion 
of free masses, it is convenient to start the analysis with the similar problem in 
electrodynamics. Let us consider the motion of a charged particle in the field of a 
monochromatic electromagnetic wave. It is known [S] that a charged particle in the 
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electromagnetic field is subject to the electromagnetic Lorentz force F given by 

F eE+ - (v X H) . (1) 

c 

The first term in Eq. (QJ is the electric contribution to the force, while the second 
term is the magnetic contribution. The ratio of the second term to the first one is of 
the order of v/c. This means that in the field of a weak electromagnetic wave, i.e. in 
the field of a wave that gives rise to a small velocity w/c <C 1 of the charged particle, 
the magnetic contribution to the force F is also small. To find the trajectory of the 
particle, one has to solve the equations 

(Px „ e / dx , \ 
"^3:t = ^e + - -xH . (2) 



dt^ ' c \dt 
The character of the trajectory depends on the wave's polarization. 

Let us start with a linearly polarized monochromatic wave of angular frequency 
uj, propagating in the direction, with the electric field E oscillating along the x^ 
axis. The trajectory of the particle, which is on average at rest, is given by |2 

^1 = — ^cosujt, a;2=0, = sin 2wt. (3) 

The particle moves in the {x^ , x^) plane along a symmetric curve, shaped as figure 8, 
with its main axis oriented in the x^ direction, see Fig. The electric component 
of the Lorentz force is dominant, and it drives the particle in the plane of the wave- 
front, while the magnetic component is responsible for the movement of the particle 
back and forth in the x^ direction. For a weak electromagnetic wave, the size of the 
a;^-amplitude is small in comparison with the wavelength: ^ ~ <C 1, but the size 
of the x'^-amplitude is even smaller: x^ ~ x^ (^) ' ratio of x^ and x^ amplitudes 
is of the order of ~ ^, where v is the characteristic value of the particle's velocity. 

A circularly polarized wave provides an exceptional case wherein the magnetic 
component of the force is inactive; the particle does not have an x'^-component of 
motion. The trajectory of the particle, which is on average at rest, is given by [3], 

x^ = — ^cosut, x"^ = --^ sin Lut, x^ = 0. (4) 

The particle moves in the {x^,x'^) plane along a circle, as shown in Fig. QJa. In the 
general case of an elliptically-polarized wave, the x'^-component of motion is always 
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Figure 1. The figure a) on the left shows the trajectory of a charged particle 
in the field of a linearly polarized electromagnetic wave, whereas the figure b) on 
the right shows the trajectory in the field of a circularly polarized wave. 

present. The projection of the trajectory onto the {x^,x^) plane forms an elUpse. The 
ellipse degenerates into a straight line or a circle for a linearly-polarized or a circularly- 
polarized wave, respectively. We shall sec below that there exists an analogy between 
the considered motion of charged particles in the field of an electromagnetic wave and 
the motion of test masses in the field of a gravitational wave. 

3. Motion of a free test mass in the field of a weak gravitational wave 

Weak gravitational waves belong to the class of weak gravitational fields: 

ds^ = c^dt^ - [Sij + hij] dx'dx^ . (5) 
A general expression for a plane wave incoming from the positive direction is given 

by 

12 

hij =Pija+Pijb, (6) 

where 

a = h+sm{k{x° +x^) + il)+), b = hxsm{k{x° + x^) + ipx), (7) 
and k = ^ = = ^, x^ = ct; V'+ and Vx are arbitrary phases. The polarization 

s 

tensors Pij {s = 1, 2) are given by 

Pij = iriimj — HjUi, Pij = —miUj — nijUi, (8) 
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where the unit vectors and rij, together with the unit vector k^/k pointing out 
in the direction of the wave propagation, form a triplet of mutually orthogonal unit 
vectors. 

One still has the freedom of turning the coordinate system in the {x^,x'^) plane 
by some angle. This transformation mixes the components of the matrix Using 
this freedom (see [Appendix A) ) one can simplify the functions Q in such a way that 
in the new coordinate system they take the form 

a = /i+sin(fc(x" +-0), 6 = /ix cos (fc(x° + x^) + V') , (9) 

where hy are arbitrary amplitudes and ■0 is an arbitrary phase. The unit 
vectors m^, Ui have the components = 1,0,0 and — 0,1,0. We shall call 
this special coordinate system a frame based on principal axes. Two independent 
linear polarization states are defined by the conditions = or /ix =0, and two 
independent states of circular polarization are defined by /i+ — hy or = —hx- 
Specifically, we will be working with the metric 

ds^ = c^dt^ - (1 + a)dx^^ - (1 - a)dx^^ + 2b dx^dx^ - dx^"^ , (10) 

where the functions a and b are given by Eq. 

The metric (jlUf) is written in synchronous coordinates. Therefore, the world 
lines ~ const are time-like geodesies, they represent the histories of free 

test masses. To be as close as possible to the framework of laboratory physics, and 
specifically to electrodynamic examples considered in the previous section, we have to 
introduce a local inertial coordinate system (t, ^). Let us associate it with the 

central world line x^ — x'^ = x'^ = 0. By definition, in the local inertial frame, and 
along the central geodesic line, the metric tensor takes on the Minkowski values and 
all first derivatives of the metric tensor vanish. The transformed metric takes on the 
form 

Cx^'^ + x^'^ + x^'^ \ 
h — 1 . (11) 

A local inertial frame realizes, as good as we can do in the presence of the gravitational 
field, a rigid freely falling "box" with a clock in it JU]. The required coordinate 
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transformation = (x") is given by [5]: 
iO = xO + ia (xi'-x2') 

(12) 

x^ = x^ - |;a {x^^ - X^^^ + ij) x^x^. 

In the above transformation, the functions a, b and their time derivatives = 

6 = are evaluated afong the world line x^ = x^ = x^ = 0. The linear and 

quadratic terms, as powers of x", are unambiguously determined by the conditions 

of local inertial frame, but the cubic and higher-order corrections are not determined 

by these conditions. In principle, transformations (|12|l can be used for all values of 

x", but they are physically useful when the values of x" are sufficiently small, that is, 

when the cubic and higher-order terms can be neglected. 

Let us consider a free test mass riding on a time-like geodesic (x^ = h, x^ = h, 

x^ = I3). Equations p2|) define the behaviour of this mass with respect to the 

introduced local inertial frame. Concretely, we have 
x^{t) = /i + i [h+li sin {tut + ip) — hxh cos {ujt + ip)] 

+ \k [h+l^h cos {ujt + V') + hy^l^h sin [ujt + ip)] , 

x^{t) = l2-\ [h+h sin(wt + V) + hy,li cos (cjt + tp)] (13) 

— i/c [h+l3l2 cos {uut -f V) — hxhh sin {uit + tp)] , 

x^(t) = ^3 - ifc [h+ - ^2^) cos {ojt + ip) + 2/ix?i?2 sin {ust + ip)] . 
Obviously, the unperturbed {i.e. in the absence of the gravitational wave) position of 

the nearby mass is (/i, hjh)- The action of the wave drives the mass in an oscillatory 

fashion around the unperturbed position. In general, all three components of motion 

are present. 

If one neglects in Eq. the terms with k, thus effectively sending the 

wavelength A to infinity, one arrives at the often cited statement about the particle's 
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Figure 2. Motion of free test masses in tlie field of a linearly polarized 
{h+ ^ 0,hx = 0) g.w. in the lowest approximation, that is, when the wave 
number k is effectively set to zero. 



motion: 



x^{t) = h + ^ [h+h sin {ujt + ip) — hxh cos {cut + 
x^{t) — I2 — ^ [h+h sin {ujt + ?/;) + hxh cos {cut + 



(14) 



Clearly, this is the analog of the electric component of motion in electrodynamics; 
moving particles remain in the plane of the wave- front. The oscillations of individual 
particles, for the particular case of a linearly polarized g.w., are shown in Fig.|21 

A circular disk, consisting of particles that are located in the plane of the wave- 
front, stretches and squeezes in an oscillatory fashion, as shown in Fig|31 

Now, let us take into account the terms with k in Eq. (|13|l . If /a = 0, these 
terms do not change the x^, x^ components of motion, but nevertheless they introduce 
oscillations along the x? direction. In analogy with electrodynamics, it is reasonable 
to call these terms the "magnetic" components of motion. The trajectories of free 
masses are, in general, ellipses, and they are not confined to the plane of the wave- 
front. In Fig. 01 (taken from 0) we show the trajectories of some individual particles. 
The "magnetic" contribution is smaller than the "electric" one. In analogy with 
electrodynamics, the major axis of the individual ellipse is small in comparison with 



Components of the gravitational force in the field of a gravitational wave 9 



Figure 3. Deformation of a disk of free test particles in the field of a linearly 
polarized {h+ 0,hx =0) g.w. in the limit of fc = 0. The two figures show the 
displacements at the moments of time separated by a half period. 




Figure 4. Solid lines show 3-dimensional motion of particles, while dashed lines 
show projections of the trajectories onto the x^,x^ plane. In the left picture the 
wave is linearly polarized (h+ ^ O.hx =0), whereas in the right picture the wave 
is circularly polarized = ft x ) ■ 

A and I, but the size of the x'^-amplitude is even smaller than the x^-amplitude. Their 
ratio is typically of the order of kl = 27r(Z/A), where I — \J Ix' + l-i' + 1-^ is the 
mean (unpertmbed) distance of the test particle from the origin. In analogy with 
electrodynamics, the "magnetic" component of motion will be present even if the 
particle is initially at rest with respect to the local inertial frame. 

A circular disk, consisting of free test particles, also behaves differently, as 
compared with the lowest order "electric" approximation. The disk does not remain 
flat while being stretched and squeezed. In addition to being stretched and squeezed it 
bends, in an oscillatory manner, forward and backward in the direction. We show 
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Figure 5. The figure shows the deformations of a circular disk of free test 
particles under the action of a linearly polarized g.w. ^ 0,/ix = 0). The 

"magnetic" contribution is responsible for the displacements along the x"^ axis. 
The two pictures show the configurations at the moments of time separated by a 
half of period. 

these complicated deformations in Fig. The level of the "magnetic" contribution is 
determined by the assumption that IjX = 0.1. This figure should be compared with 
Fig. El 

4. Equations of motion from the geodesic deviation equations. 

In this section we shall consider the motion of free masses from a different perspective. 
The geodesic deviation equations will be used in order to derive the analog of the 
Lorentz force Eq. (^) and the analog of the Newtonian equations of motion Eq. Q . 
We shall see that the "magnetic" component of motion is contained in the higher-order 
geodesic deviation equations. 

J^.l. Geodesic deviation equations in general 

The derivation of the geodesic deviation equations is usually based on a 2-parameter 
family of time- like geodesies x^{t, r). For each value of r, the line x^{t, r) is a time- like 
geodesic with a proper time parameter r. The vector is the unit tangent vector to 
the geodesic, and is the "separation" vector between the geodesies (see Fig. 01: 

dx^ Ox^ 

u'^{T,r) = ——\r=const , n'^ {t, r) = ——\r=const ■ (15) 

OT or 
Let the central (fiducial) geodesic line correspond to r = 0, and a nearby geodesic 
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to r — rQ. For small tq, one has: 

x^{r, ro) = x^{r, Q) + r,^ + -r^^ + 0(rg). (16) 

Then, in the lowest approximation in terms of tq, the geodesic deviation equations are 
given by pil] : 

^ = K,,n-u^n-', (17) 

where i?^^^ is the curvature tensor calculated along the geodesic line r = 0, and ^ 
is the covariant derivative calculated along that line. 

To discuss the "magnetic" component of motion in the field of a gravitational wave 
we will need the geodesic deviation equations extended to the next approximation. 
These equations were derived by Bazanski [111 . A modified derivation can be found 
in jl2) . In the required approximation, one needs the information on the second 
derivatives of x'': . These quantities do not form a vector. It is convenient to 
introduce a closely related vector w^: 

-'' = ^=->" = ^+r>"n^. (18) 
This vector obeys the equations JT] (see also \V2\]: 

^ - K,,^"^'^' + {k,-,s ~ + 4i?;:,,-^^n^. (19) 

To combine equations (|17() and H19|l into a single formula, valid up to and including 
the terms of the order of Tq, it is convenient to construct a vector N^: 

= r^n^ ^\^rlw^ . (20) 

Taking the sum of Eq. (|17() . multiplied by ro, and Eq. I|19|l multiplied by ro/2, we 

obtain the equation whose solution we will need, 
rP-N^ 1 / \ 

ap7 2 \ ctP7;o 'y6oi.;lj j 



+ '^Kpy^N^ + 0{rl). (21) 

In terms of , the expansion (|16|) takes the form 

x^'{T, ro) = x^{t, Q) + - T^'^^N'^N^ + 0{rl). (22) 

Formula H22I) shows that in the local incrtial frame (in which case F^^ = along the 
central geodesic line) the spatial components TV' will directly give the time-dependent 
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Figure 6. Deviation of two nearby geodesies in a gravitational field. Line (a) 
represents the central geodesic r = and line (b) represents the nearby geodesic 
r = ro- 11^ is the unit tangent vector to the central world line, n^* is the tangent 
vector to the curve t = const. The lines (d) and (c) are the corresponding first 
and second order approximations to the exact nearby geodesic (b). 



positions of the nearby particle. According to Eq. (|21(l . these positions include the 
next order corrections, as compared with solutions to Eq. I|17() . 

In Fig we show the successive approximations to the exact position of the 
nearby geodesic (6). The line (d) represents the first order approximation (i.e. 
x'^{T,ro) = x'^(t, 0) + ro^^ ), the line (c) takes into account the second order 
approximation {i.e. x'^(t, rp) — x''(t, 0) + ?'o^|r + i^o^^ )' according to ((TB|l . 

4-2. Geodesic deviation in the field of a gravitational wave 

We now specialize to the g.w. metric (0), ©, where a and b are given by ©. We 
take account only of linear perturbations in terms of g.w. amplitude h. The first 
particle is described by the central time-like geodesic x^{t) = 0, its tangent vector is 
It" = (1, 0, 0, 0). The second particle resides at the unperturbed position x*(0) = P and 
has zero unperturbed velocity. It is assumed that the local inertial frame is realized 
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along the central geodesic. The objective is to find the trajectory of the second particle 
using the geodesic deviation equation H21|l . 
The deviation vector N"^ has the form 

N'{t)^l' + Cit), (23) 

where ^*(t) is the perturbation caused by the gravitational wave. The choice of the 
local inertial frame allows one to replace all covariant derivatives in Eq. (|21|) by 
ordinary derivatives. In particular, covariant time derivative is being replaced by 

,2 

ordinary time derivative ^3^- In the lowest approximation, Eqs. H21|l reduce to Eqs. 
IjlTfl and specialize to 

I.,,- 92 1 



1 (9^ 
2' dt' ^ 



1 . 2 . 

P ]h+ sin(wi + Tp)+P)hx cos(wt + ip) 



(24) 



dt^ 2 dt^ ^ 2- 
As expected, the relevant solution to this equation coincides exactly with the usual 
"electric" part of the motion, which is given by Eq. (|14|l . 

In order to identify the "magnetic" part of the gravitational force, one has to 
consider all the terms in Eq. (|21|l . Since DN" /dr is of the order of h, the third term 
in Eq. H21|l is of the order of and should be neglected. Working out the derivatives 
of the curvature tensor and substituting them into H21I) . we arrive at the accurate 
equations of motion: 

d^C 



2 dt^ ^ ^' 2 dt^ \dx^ ^ 2 dx3 ^ 7 



(25) 



The second term of this formula is responsible for the "magnetic" component of motion 
and can be interpreted as the gravitational analog of the magnetic part of the Lorentz 
force (HJ. 

Specifically, in the field of the gravitational wave IjlUf) , the full equations of motion 
H25|l take the form: 



rf'e 1 2n 

-d^^r' 



P ih+ s\\\{uot + V^) + P j/ix cos(wt + tp) 



2 



P kjh+ cos{iLjt + ip) — P kjhx sin{ujt + -0) 



This equation clearly exhibits two contributions: 



pi 



(m)- 



(26) 



(27) 
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The "electric" component of the gravitational force F^^-^ /m is given by the first term in 
Eq. (jSnj. The second term - the "magnetic" component F^j^^/m of the gravitational 
force - can be written, demonstrating certain analogy with electrodynamics, in the 
form involving the (lowest order non-vanishing) velocity ^ of the test particle: 



pi 



The right hand side of Eq. (|26|l is exactly the acceleration ^jp- which can be 
derived by taking the time derivatives of Eq. I|13|l . Not surprisingly, by integrating 
equations of motion (|26|) . one arrives exactly at the time-dependent positions of 
the particles, Eq. H13|) . which we have already derived by the direct coordinate 
transformation. Therefore, the gravitational Lorentz force, identified above, leads 
exactly to the expected result. 

It should be noted that Eqs. H26|) depend only on the so-called TT-components 
of the g.w. field. This happens not only because we have explicitly started from 
them in Eq. (|10|l . Even if we have started from the general form of the g.w. field, 
which includes also the non-TT components, we would have ended up with equations 
containing only the TT-components. This happens because Eqs. (|21|) involve the 
curvature tensor (and its derivatives) in which the non-TT components automatically 
cancel out. 

It is interesting to compare the components of the gravitational force derived 

here with what would follow from the concept of "gravitomagnetism" . In general, 

the concept of "gravitomagnetism" is a helpful analogy which was successfully used 

in studies of stationary gravitational fields ^31) EH- However, its application to the 

gravitational-wave problem considered here requires certain care. For example, in 

the local inertial frame, the leading terms of the equations of motion derived in the 

framework of "gravitomagnetism" |14 | read (in notations consistent with this paper): 
<PN^ dN^ , 

-^-V + ^Rm^N'^- (29) 

This equation should be compared with our Eq. H21|) . also specialized to the local 
inertial frame: 

d^ 1 dN-^ 
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Clearly, the last term in both equations is common, and it resembles the magnetic 
part of the electromagnetic Lorentz force. However, as was shown above, for particles 
which do not have large unperturbed velocities and are (on average) at rest in the 
local incrtial frame, this term is quadratic in h and should be neglected. At the same 
time, the second term in Eq. H3U|) (which we provisionally call "magnetic") is linear 
in h and cannot be neglected. Regardless of terminology, the correct results exhibited 
in Eq. (|13|l can only be obtained if one proceeds with Eq. H30() and not with Eq. 129() . 

5. Variation of the distance between test masses 

We have used a local inertial frame and completely specified the time-dependent 
positions of test particles acted upon by a gravitational wave. As was already 
mentioned, this description is as close as possible to the description of laboratory 
physics. Having answered all the questions with regard to particles's positions, we 
can now discuss the variation of distances between them. We are mostly interested 
in the distance between the central particle and the particle located, on average, at 
some position {h,l2,l3)- This is a model for the central mirror and the end-mirror 
placed in one of the arms of an interferometer. We will later use these results for the 
derivation of the response function of the laser interferometer. 

In the local inertial frame, metric tensor (|ll|l has the Minkowski values up to 
small terms of order of h{j)^. The Euclidian expression 

d{t) = \/5"%^?2~^ ^ Q {hl{l/\f) , (31) 

gives the distance between particles, which is accurate up to terms of the order of 
hi and hlj inclusive. Obviously, we neglect the terms quadratic in h. Denoting 
— li + Ax', one gets 

d{t) ^l + j ('lAx^ + hAx^ + kAx^) . (32) 

Using the time-dependent positions (|13|l . we obtain the distance d{t) with the required 
accuracy: 

d{t) = 1 + [h+ih'^ - h^) sin {^^t + V') - 2hxlil2 cos {ujt + ^p)] 
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+ [h+{li'^ - h^) cos {ut + -0) + 2h^lil2 sin {ut + -0)] + O {hl{l/Xf) . 

(33) 

Clearly, the first correction to / is due to the "electric" contribution, whereas the 
second correction to / is due to the "magnetic" contribution. 

According to Eqs. (|13|l . the "magnetic" component of motion is present even if 
the mean position of the second mass is such that Z3 — 0. However, this motion is 
in the direction orthogonal to the line joining the masses and therefore it does not 
lead to a (first order in terms of h) change of distance between them. This fact is 
reflected in Eq. (|33|l in the form of disappearance of the "magnetic" contribution to 
the distance when ^3 = 0. In other words, "magnetic" contribution to the distance is 
present only if the interferometer's arm is not orthogonal to the wave's propagation. 

The distance (|33|l was calculated in the local inertial frame. It was assumed that 
l/\ ^ 1. It is important to show that the approximate expression (|33ll follows also 
from exact definitions of distance. One of them is based on the measurement of time 
that it takes for a light ray to travel from one free particle to another and back. [This 
is a part of a more general problem of finding light-like geodesies in the presence of a 
weak gravitational wave 2ni>Cni> EI-] This definition is applicable regardless of the 
relationship between I and A and does not require the introduction of a local inertial 
frame. If a photon is sent out from the first particle-mirror at the moment of time 
to, gets reflected off the second particle-mirror, and then returns back to the first 
particle-mirror at t2, the proper distance d{t) between the mirrors at time t is defined 
as 

d{t)=c'-^, (34) 

where t = (to + i2)/2 is the mean time between the departure of the photon and its 
arrival back. 

In the field of the gravitational wave H1U|I . the light rays, ds^ = 0, are described 
by the equation 

c^dt^ = (1 + a)dx^^ + (1 - a)dx^^ - 2b dx^dx^ + dx^^ . (35) 
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Let the (unperturbed) outgoing light ray be parameterized as 



X 



° — rfo + It, = liT, x^ = I2T, x^ = Z3T, (36) 



where the parameter r changes from to 1. Then, according to Eq. (|35|l . we have 
along the ray: 



cdt — I 



1 + ^«W^-Kr)^ dr. (37) 



Integrating both parts of this equation, we can find the time t of arrival of the photon 
to the second particle. The calculated time includes the g.w. corrections proportional 
to /i+ and hx- Similarly, the (unperturbed) reflected light ray can be parameterized 
as 

x^ ~ ct + It, x^ = li — liT, x^ ~ I2 — I2T, x^ ~ I3 — I3T, (38) 

where r is again changing from to 1. A similar integration of Eq. H37II allows us to 
find the time t2 of arrival of the photon back to the first particle, including the g.w. 
corrections. 

Combining two pieces of the light travel time, we derive the exact formula for the 
distance: 



dit)^l + l-^ 



72 _ ;2 
'1 '2 



$c + hxhh'^s 



(39) 



where 



= — TT, TT [cos (wt + + ^^3) — COS (cjt + " ^0] 

K(l + h) 

+ , ,, ^ , , [cos (wt + + kh) - cos (ujt + V + kl)\ , (40) 
k{l - h) 

and $s is obtained from $c by the replacement of all cos-functions with sin-functions 
of the same arguments. Exactly the same formula (|39() follows also from the direct 
joining of the perturbed light-like geodesic lines, derived in reference jl5| . 

Formula H33|l is sufficient for ground-based interferometers, for which the 
condition Z <C A is usually satisfied. Formula (|39|l is appropriate for space-based 
interferometers for which the above condition is not satisfied in the higher-frequency 
portion of the sensitivity band (see, for example, J7], JS], ^U). However, it 
is important that formula l|33f) . including its "magnetic" terms, follows also from 
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exact definition (|39|l . when the appropriate approximation is taken. Assuming that 
k{l + I3) <ti 1 and k{l — Z3) ^ 1 and retaining only the first two terms in the expansion 
of $c and $s, one derives Eq. IpHjl from Eq. As expected, the "magnetic" 

contribution to the distance is a universal phenomenon. It is most easily identified 
and interpreted in the local inertial frame, but conclusions about the distance do not 
depend on the introduction of this frame. 

It is known that there is no unique definition of spatial distance in curved space- 
time. We have considered the definition based on measuring the round-trip time of a 
light ray. One more definition is based on measuring the length of a spatial geodesic 
line joining the particles at a fixed moment of time. It can be shown that this definition 
leads to the formula 



p - p ~ 



(41) 



where 



2 

$c = - -n- [cos (LLjt + tp + kh) - cos iujt + V)] , (42) 
fc/3 

and $s is obtained from $c by the replacement of all cos-functions with sin-functions 
of the same arguments. In general, the distance (|41() differs from the distance H39() . 
However, they do coincide in the first order approximation in terms of small parameter 
kl. This can be shown by expanding Eq. I|41|l in powers of kl and retaining the linear 
terms. It is satisfying that the exact definitions H39|) and \W\ lead precisely to the 
Euclidean result %6'6\ in the appropriate approximation. 

6. Response of an interferometer to the incoming plane gravitational wave 

Laser interferometer measures the difference of distances travelled by light in two arms. 
We describe interferometer in the local inertial frame, with the origin of the frame 
at the corner mirror. The unperturbed coordinates of the end-mirrors are given by 
{lf\ if^, if^), where a = 1, 2 labels the arms. We consider an interferometer whose 
unperturbed arms have equal lengths, I = \l l\ ' +I2 + ^3 : and the arms are 
orthogonal to each other. A detailed introduction to interferometers in gravitational- 
wave research can be found in 1201) P)- 
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It was shown in the previous section that the distance variation in one of the arms 
is given by Eq. 1)33(1 . Then, the response of a 2-arm interferometer is given by 

Ad(i) = d(t)(i) -d(t)(2). (43) 

The derivation of formula (|33|l is based on the coordinate system adjusted to the 
gravitational wave. Specifically, it is assumed that the axis is the direction from 
which the incoming plane wave propagates, while the x^,x'^ directions are defining 
the principal axes of the wave, see Eq. H1U|) . Then, the response function H43|) is 
characterized by six free parameters {l[^\ 4^'i 4^')- Among these six parameters 
only four are independent (one of which is I), because the arms have equal unperturbed 
lengths I and arc orthogonal to each other. 

When it comes to the observer, it is more convenient to associate coordinate 
system with the interferometer's arms, rather than with one particular wave. Let the 
observer's coordinate system {X^,X'^, X^) be chosen in such a way that the arms are 
located along the X^ and X^ directions. Then, the response function is characterized 
by I and three angles : 8, $ and 5*. The angles 8,$ describe the direction of 
a particular incoming wave, and the third angle 5" describes the orientation of the 
principal axes of the wave with respect to the observer's meridian (see, for example, 
12], [201 J 123 )• The relationship between these two coordinate systems is shown in 
Fig. [7| and is described in detail in [Appendix B| 

The parameters (4^\ 4^"*' 4^'') are expressible in terms of the parameters 8, $, ^ 

according to the relationships (see [Appendix Bj ): 

= ;(cos$cos8cos* - sin $ sin ^P), 
4^^ = ;(- sin^cos* - cos $ cos 8 sin*), (44) 
4^^ = ;(cos$sin8), 

and 

4^^ = ;(sin$cos8cos* + cos$sinvIf), 

4^-' = ;(cos$cos* - sin$cos8sinvIf), (45) 
4^^ = ;(sin$sin8). 
Using H44|l and <|45|l in (|43|l . the response function can be written in the form 

Ad{t) =l[h+ {Ff (8, *) sin(tjt + i;) + {kl)F^{e, $, *) cos(wt + i/;)} 

+/ix {F^{Q, cos(wi + V) + {kl)F^{Q, $, ^) sm{ujt + V)}] , (46) 
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Figure 7. The relationship between the coordinate system {x^,x^,x^) 
adjusted to the wave and the coordinate system {X^,X^,X^) adjusted to the 
interferometer. The two arms of the interferometer lie along the axes and X^. 



where the "electric" {E) components of -F+.x (©, 'I', are given by 



Ff(e,$,*)= cos 2$ 



F^(e,$,*) = -cos2$ 



1 + cos^ e 



1 + cos^ e 



cos 2^1^ - sin 2$ cos Q sin 2^, 



sin 2* - sin 2<i> cos 8 cos 2^, 



(47) 
(48) 



and the "magnetic" (H) components of F+.x (©, ^, ^') are given by 

' 1 + cos^ e 



(9,$,^-) = isine 



cos 8 + sin 2<f> 



sin 2$ (cos $ + sin $) cos 8 sin 2"^! 



(8,*,^') = - sine 



cos 8 + sin 2<f> 



1 + cos^ 8 



(cos$ - sin<i>) cos 2^1/ 



(cos <I> — sin $) sin 2^' 



(49) 



(50) 



+ sin 2$ (cos $ + sin $) cos 8 cos 2\E' . 

Our response function is more accurate than the previously derived expressions 
| 20) . |21| . |22) . because our Eq. I|4t)|) includes the "magnetic" contribution. In equation 
H4()|l it is given by the terms proportional to the factor {kl). 

In general, the response function H46() contains two independent polarization 
amplitudes, h+ and h^- To simplify the analysis of Eq. H46|l . we will separately 
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consider circularly polarized {h^ = ±h-^), and linearly polarized (/ix = or = 0), 
waves. We start with the analysis of the response function as a function of time, 
assuming that a circularly polarized or a linearly polarized wave arrives from a fixed 
direction 8, $ on the sky. 

In the case of a circularly polarized wave {h^ = ±/ix) we obtain 



Ad(t) ^lh+ 

+ + ■ Fr + F 

where the phase shift Aijj is given by 



sin(wt + ?/'± AV'), (51) 



tan (AVj) 



F^ \ .../F^F]'-F^F"\1 



X 



1 ± (kl) 



+ ^ X ^ X 

pEpE 



(52) 



Clearly, the inclusion of "magnetic" terms changes the amplitude and the phase of 
/S.d{t). As an illustration, we show in Fig. |Slthe response of an interferometer, as a 
function of time, to a circularly polarized wave = /ix, coming from the direction 
G = 27r/3,<I> = 7r/3 (we have also set ^l* = to fix the phase of the response function; 
the amplitude of the response function docs not depend on ^f). The dashed curve 
shows the "electric" response alone, while the solid curve shows the total response, 
including the "magnetic" part. We have taken l/X — 0.1. 

The response of an interferometer to any elliptically polarized wave is qualitatively 
similar to Fig. |H1 that is, in general, the "magnetic" component contributes, both, 
to the amplitude and to the phase of Ad{t). However, "magnetic" contribution to 
the amplitude may be very small for linearly polarized waves. If one puts hx ~ or 
/i+ = in Eq. H46|l . one finds that the correction to the amplitude Ad{t) is quadratic 
in (kl)'^ (and, hence, can be neglected) with the exception of directions on the sky 
where F^ or F^ vanish. Specifically, for the linearly polarized wave hx — 0, we 
obtain 



Ad{t) = / h+ JFf + {kl)^F^ sin {cut + i/j + A^p), (53) 



where the phase shift Aip is given by 

-pH 



tan(AV') = fc/(^-^j . (54) 
(The case when — is given by the above expressions in which all the "plus" indices 
are replaced by "cross" indices, and the sin-function in expression H53|) is replaced by 
a cos-function.) 
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response 




Figure 8. A typical response of an interferometer, as a function of time, to 
the monochromatic circularly polarized gravitational wave coming from a fixed 
direction on the sky. The solid line shows the total response, while the dashed 
line is purely "electric" part. 

We now turn to the amplitude of Ad{t) as a function of the angles 8,$. It 
is usually called the beam pattern. In general, the amplitude of Ad{t) depends 
also on but we will consider, for simplicity, circularly polarized waves, in which 
case the parameter ^ does not participate in the amplitude of Ad{t). We consider 
circularly polarized waves, with one and the same amplitude hn, coming from arbitrary 
directions on the sky. The beam pattern is shown in Fig. O The left figure shows the 
purely "electric" contribution, whereas the right figure shows the total beam pattern, 
with the "magnetic" contribution included. We have taken l/\ = Q.l. It is clearly 
seen that the "magnetic" component breaks the characteristic quadrupole symmetry 
of the "electric" beam pattern (more detail - in the next Section). 

7. Response function in terms of spin-weighted spherical harmonics 

The response function H4t)|l allows an elegant invariant representation in terms of 
the spin-weighted spherical harmonics. This description is also helpful for physical 
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Figure 9. The amplitude of the interferometer's response to circularly polarized 
waves. The graphs are normalized in such a way that the amplitude is equal 1 
for = 0. The left figure ignores the "magnetic" effect, whereas the right figure 
shows the total response. 



interpretation of the response function. 

First, let us introduce the amplitudes of circularly polarized gravitational waves, 



and the complex functions 
Fn = (Ff 



{kl)F^)+^{F^ + {kl)F^), 



(55) 



(56) 



Fl = (Ff - ikl)F^) + ^ {-F^ + ikl)F^) . 
In terms of the introduced notations, the response function ()46|l can be identically 



rewritten as 



Adit)^l[-- 



[HrFr + HlFl) - e 



{HrF*^ + HlFD 



(57) 



Clearly, the Ad{t) has been separated in two components associated with the left (L) 
and right (R) polarization states of the gravitational wave. Only the L{R) component 
of the response function is present, when the R{L) g.w. amplitude is set to zero. 
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Using the explicit expressions (|17|l . and (|5n|) . one can show that the 

newly defined functions Ffj, Fl can be rearranged to read 

Fn = Fn{Q,^)e''''', = i?i(e, $)e2^*, (58) 

where Fl and Ffi are given by 



(1 + cos 9)' e^^* + (1 - cos 9)' e~^'^ 



ikl 
32 



(1 - i) (1 + cos 6) sin 9e-^'* + (1 + i) (3 cos 9-1) (cos 9 + 1) sin 9e*' 



+ {l-i) (3 cos 9 + 1) (cos 9-1) sin 9e"'* + {l + i){l- cos 9)^ sin 9e 



(1 - cos 9)' e^'* + (1 + cos 9)' e'^"^ 
ikl 



-3i<S> 



(59) 



32 



(1 - i) (1 - cos 9)'' sin 96^** + (1 + i) (3 cos 9 + 1) (cos 9-1) sin Qe' 



+ (1 - i) (3 cos 9-1) (cos 9 + 1) sin 96"'* + (1 + i) (1 + cos 9) sin 9e 



(60) 



For a given direction (9, $) on the sky, the functions Fj^ and F^ transform under 
the rotation of ^P, specified by the transformation ^ ^I* + ^'o, according to Fn 
fi{e~^**° and Fl — > i^i,e^**°. A function s/(9,$) is said to be a spin-s function, if 
under the transformation '5 ^ 'f + ^E'o it transforms like s/(9, <f>) s/(9, (f>)e^*'**". 
Therefore, our functions Fl and represent the s = — 2 and s = +2 fmictions, 
respectively. 

A scalar function on a 2-sphcre can be expanded over a set of ordinary spherical 
harmonics YJ^{Q, $), which form a complete and orthonormal basis. These harmonics 
are not appropriate for the expansion of spin-weighted functions with s ^ 0. The 
spin-weighted functions can be expanded over a set of harmonics called the spin- 
weighted spherical harmonics sYmiQ,^) ESI. EH. HHIi QHI' The spin-weighted 
spherical harmonics satisfy the completeness and othonormality conditions similar 
to those of ordinary spherical harmonics. The functions sY,l^{Q,^) can be derived 
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from 1^(0, $) according to the rules: 

''B'Y^, {0<s<l) 

1/2 
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1/2 



(61) 



^11, (0>s>-0 

where the differential operators S and 8 [acting on a spin s function s/(6, <&)] are 
given by 

g ,/(e, $) = - sin^ e + sin-^ e ,/(e, $), 



(62) 



lae ^ sine 

8 «/(e, $) = - sin-« e (^ - ^ ^) sin« e ,/(e, $). 

Using the spin- weighted spherical harmonics, the functions jFl and jFn can be 
expanded as follows: 



/ " 



^)_,y33(e,$) + ii±^_2n3(e,$) 



+ 



^/3 ■ 
(1-i) 

a/5 ■ 



a/5 

.2F_\(e,$) + ii±!l_2y33(e,$) 



(1-^ 

a/3 



V3 



2y33(e,$) + i^2lf(0,*) 



(63) 



'V_\(e,$) + ii±^2F^(e,$) 



(64) 



V5 ^^-1^^'^^ V3 
As one could expect, the "electric" components of the response function are described 
by the spin 2 quadrupole (1=2) terms ±2^, whereas the "magnetic" components are 
described by the spin 2 octupole (1=3) terms ±2^- The "magnetic" components 
can be viewed as the higher-order terms in the multipole expansion of the response 
function. 

8. Astrophysical example 



The "magnetic" component in the motion of free masses will have serious practical 
implications for the current gravitational wave observations. For example, the LIGO 
interferometers have the arm length of / — Akm and arc most sensitive to g.w. in the 
interval of frequencies between 2>QHz and IQ'^Hz. This means that the "magnetic" 
component, whose relative contribution is of the order of kl, may provide a correction 
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at the level of 5% to the response function in the frequency region of 600Hz, and up 
to 10% in the frequency region of l200Hz. This contribution may significantly affect 
the determination of the source's parameters. 

The high frequency region, QOQHz— l200Hz, will be extensively studied with the 
help of the "narrow-band tuning" of advanced interferometers [S] . The signal recycling 
technique allows reshaping of the sensitivity curve in such a way that it reduces noise 
within a chosen narrow band, while weakening sensitivity in other frequency regions. 
This region of high frequencies is populated by various periodic and quasi-periodic 
astrophysical sources, such as neutron stars in low-mass X-ray binaries, including Sco 
X-1, slightly deformed rotating neutron stars, neutron star and black hole binaries in 
the last moments of their inspiral. For a more detailed list of high frequency sources 
and the prospects of their detection see |S]. 

To better understand the role of the "magnetic" components in estimation of the 
g.w. source parameters, we shall briefly consider an idealized example of a compact 
binary system. Let the binary consist of objects of equal mass M orbiting each other 
in a circular orbit of size 2r. The distance of the binary to the observer is R. For 
simplicity, we assume that the orbital plane is orthogonal to the line of sight defined 
by 8;$. Then, the gravitational wave at Earth has the form of Eq. ^ with the 
amplitudes given by [20], 0]: 

h+ = h^= Hr = ^Mr^cj^, (65) 

where the g.w. angular frequency twice the orbital frequency. The angle 

\1/ can be taken as 5* = 0, and the phase ip is the angle between the observer's meridian 
and the line joining the components of the binary at some initial moment of time. 
The response of the interferometer to the incoming wave is determined by Eq. 

Ad{t) = Ihn |cos2$ ( ^^'^2^ ^ ^ sin 2$ (cos$ + sin*) cos9sin9 



X sin {utt + ip) 

sin2(|>cos9 - ^— ( cos^ Q + sin 2$ ( " ' ) ) (cos$ - sin<i>) sin 6 
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X cos {ujt + ip) 



(66) 



The terms with ojI represent "magnetic" contribution. If the observed Ad{t) is 
incorrectly interpreted as the "electric" contribution only, then the parameters of 
source (for example, its mass 2M) would be estimated from the relationship: 



Clearly, this would have resulted in a significant error in the estimated M (of the order 



response function, consisting of "electric" and "magnetic" contributions. A concrete 
value of the error depends on the direction 0, 

9. Conclusions 



In this paper we have considered the motion of free test particles in the field of a 
gravitational wave. We have shown that this motion is similar to the motion of charged 
test particles in the field of an electromagnetic wave. Using different methods we have 
demonstrated the presence and importance of what we call the "magnetic" components 
of motion. Regardless of interpretation and terminology, these terms contribute to the 
variation of distance between the interferometer's mirrors and, hence, they contribute 
to the total response function of the interferometer. The "magnetic" contribution 
must be taken into account in advanced data analysis programs. 
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Appendix A. The principal axes for a monochromatic gravitational wave 



In general, the metric components are given by Eq. jT)), where ip+ and ipx are arbitrary 
independent constants. However, it is always possible to do a rotation in the {x^ , x'^) 
plane, such that in the new coordinate system x^,x'^ the difference between ipx and 
^-(- is exactly ^. The required transformation 



= x-^ cosC + a:^ sin^, 



= — a;-^ sin ^ + x^ cos 



has the angle C such that 



C — — arctan 
^ 2 



hi - hi 



± 



(A.l) 



,(A.2) 



2/i+/ix cos (?/>+ — psix) J y \2h-^-hx cos (1^+ — ipx) , 
The rotation angle ^ is well defined, except for a circularly polarized wave, in which 
case the difference between tpx and ip^ is always ^. In this new coordinate system 
(5^,5;^) the metric components take the form 

a = h^ sin (fc(2;" + x^) +4'), 

where the amplitudes /i+ and hx are given by 



tlx cos {k{x" + x^)+il;), (A.3) 



— (/i^ cos^ 2( + 2h+hx sin 2^ cos 2^ cos 



hy^ = (h\ cos^ 2C - 2h+hx sin 2^ cos 2C cos - ipx ) + h\ sin^ 2C) ^ , 
and the phase ip is given by 

/i-f sin'!/'_|_ + hx sin-f/^x tan2C 



ip = arctan 



/i+ cos ip-i- + hx cos 'i/'x tan 2^ 



(A.4) 



(A.5) 



Components of the gravitational force in the field of a gravitational wave 29 
We call the coordinates and x'^, in which the metric components take the form HA.3() . 
the principal axes of the gravitational wave. In the text, we assume that the above 
transformation has already been performed. So, without loss of physical generality, 
we put Tp^ — ~ § = V'l E^nd we do not use the tilde- symbol. 

Appendix B. Coordinate transformation 

A general gravitational wave H1(J|) , @ is described in the coordinate system (x^ , , x'^) . 
For the observer, it is convenient to use the coordinate system (X^, X^, X^) in which 
the arms of interferometer are located along the and X'^ axes. In this coordinate 
system, the incoming wave is characterized by the direction O, $ and the angle 4" 
between the observer's meridian and the wave's principal axis, as depicted in Fig. d 
The relationship between these two coordinate systems is given by the transformation 

(B.l) 

where the orthogonal rotation matrix R is 
f cos $ cos O cos ^ — sin $ sin ^ — sin <1> cos ^ — cos $ cos O sin 'I' cos "1> sin B \ 

sin $ cos 6 cos * + cos $ sin * cos $ cos 'I' — sin $ cos 6 sin 'J' sin $ sin 6 . (B.2) 

\ — sin O cos ^ sin 9 sin ^ cos O / 

In the observer's frame, the end-mirrors have the coordinates [X^ = l,X^ = 

0,X^ = 0) and {X^ =0,^2 = l,X^ = 0). Therefore, in the g.w. frame, the 

coordinates of the end-mirrors are 

/ cos $ cos 6 cos 4* — sin $ sin 'i> \ 



X^ \ 




' x^ 


x^ 




x^ 


X -' J 




, x^ 



(1) - R-1 



R-^ I 



and 



if^ = I 




= l 




= l 



— sin $ cos 5* — cos $ cos 6 sin 
cos $ sin 

/ sin <i> cos 9 cos + cos $ sin * ^ 
cos $ cos ^' — sin $ cos 9 sin 5' 



(B.3) 



(B.4) 



y sin $ sin 9 / 

These formulae allow us to parameterize the positions of the end-mirrors in terms of 

l and the angles 9, $, ^I'. 



